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We estimate the contribution of the long-wavelength el-TO-ph interaction and discuss the effect it 
has on electron pairing in doped polar systems like SrTi03 and PbTe. The analytical and numerical 
results presented in the study indicate that the el-TO-ph interaction tends to contribute little to 
the total strength of electron-phonon coupling in these and related materials. To explain this fact 
we consider possible reasons why the effect of the polar long-wavelength transverse optical phonons 
on the coupling constant A is far less than one might suppose. 
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I. INTRODUCTION 

A classical representative of polar crystals SrTiOs, 
which is the prototypical among insulating ABO3 per- 
ovskite oxides, is a well-known material with exceptional 
dielectric properties that made its attractive for a num- 
ber of electronic applications. The ion-covalent bonds, 
the long-range dipole forces and extremely close prox- 
imity to the ferroelectric state determine SrTiC>3 as a 
model system with a great potential for scientific inves- 
tigations. The electron-doped SrTiC>3 represents an ex- 
ample of a low carrier density system with a metallic-like 
conductivity. Doping into insulating state of SrTi03 is 
performed by electron donating ions like La 3+ for Sr 2+ or 
Nb 5+ for Ti 4+ which transfer the doped electrons into the 
originally empty d° configuration 1]. Experimentally ob- 
servable superconductivity in the electron-doped SrTiC>3 
has been assumed to be described by the plasmon-polar 
optic phonon mechanism 

The polar-coupling electron-phonon system in SrTiOs 
has attracted a substantial amount of interest 0-Q ■ The 
most important reason for this is that in such materials 
there exists a strong coupling between electronic and lat- 
tice degrees of freedom. Correspondingly, the properties 
of conductivity electrons are strongly affected by local 
lattice properties. The main motivation for the present 
study is related to previous work, in which we showed 
that in polar crystals like SrTiOs the strength of the in- 
teractions between valence electrons and transverse opti- 
cal (TO) phonon modes is especially large @, [1]. Appli- 
cation of these results to the case of doping will allow us 
to have a better understanding of how electron dynam- 
ics is affected by the strong el-TO-ph interaction. Note 
also that a good knowledge about the roles of the coex- 
isting interactions responsible for various physical prop- 
erties is valuable to further advances in polar and related 
systems. Therefore, it would be desirable to examine 
in detail whether and under which conditions does the 
el-TO-ph interaction contribute to the total electron- 



phonon coupling constant. The accomplishment of this 
task is the scope of the present communication. 



II. MODEL 

Consider a linearly coupled system of electrons and 
polar (infrared-active) TO lattice vibrations. Assume 
that the band structure is represented by the set of one- 
particle Bloch wave functions | crk> and energies E a (k) 
which can be determined within the framework of DFT- 
LDA methods with the correct accounting for the quasi- 
particle effects. The relevant Hamiltonian, which de- 
scribes the dynamic mixing of electronic states of the 
valence (tr = 1) and conduction band (a = 2) caused by 
the TO phonons, reads 



H, 



el—ph 



J2daa 

T,a' , j k,q 
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where the quantities g aa i (qj) describe the el-TO-ph in- 
teraction, a + (a) are the creation (annihilation) opera- 
tors for the states | crk> \ a=1 2 , N is the number of unit 
cells; the TO phonon modes are characterized by the nor- 
mal coordinate Uqj , the wave vector q and the vibration 
branch j . 

The consistent microscopic analysis of the model (fTJ) 
is presented in [8[. Some relevant characteristics of this 
model are summarized as follows: (1) Eq. (JTJ) is mostly 
related to polar crystals in which the difference in long- 
range fields given by longitudinal and transverse optical 
phonon modes is essential 0. (2) For q — 0, the inter- 
action of band electrons and TO lattice vibrations has 
completely the interband nature (the intraband (cr = a ) 
terms in Eq. (fTJ) vanish due to inversion symmetry [H,|8j]). 
(3) The model represents the Frohlich-type interaction 
which may give rise to an asymmetric charge distribu- 
tion Q. (4) At the macroscopic level, the strength of 
the el-TO-ph interaction is the direct effect of coupling 
between the dipole polarization associated with the po- 
lar TO mode [10| and the resulting shift of the electron 
density Q. Note also that the electronic and structural 
significance of the el-TO-ph interaction is incorporated 
within the vibronic model (e.g., Q). An example here 
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can be a family of the ferroelectric ABO3 perovskite ox- 
ides, in which hybridization (|TJ) involves strong dynamic 
mixing between the O 2p and the d° (Ti 4+ , Nb 5+ /Ta 5+ , 
etc.) states caused by the TO Fi u soft vibrations [lTI.[T^|. 
Estimates of the el-TO-ph interaction constants for a 
wide number of polar compounds are reported in 0] . 



III. THEORETICAL ANALYSIS IN THE 
LONG- WAVELENGTH LIMIT 

A widely-used approach for estimating the electron- 
phonon coupling strength is formulated in terms of 
the quasiparticle mass enhancement, m*/m, due to the 
electron-phonon renormalization factor 1 + A (e.g., (13 , 
0): 



ui^j and the additional term Aw^- 



which accounts for 



m* jm = 1 + A . 



(2) 



Here A is a dimensionless effective constant which charac- 
terizes the total strength of electron-phonon coupling. As 
well-known [TU, [lij , transition temperatures in conven- 
tional superconductors depend directly on A . According 
to common definitions (the representation of A as a sum 
over mode coupling constants) one can define the partial 
coupling constant Ato which is related to the correspond- 
ing contribution of the el-TO-ph interaction to A . 

Being based on the model ([T]), our goal is to deter- 
mine which characteristic values of Ato (J) f° r the Jth TO 
mode correspond to the el-TO-ph interactions (which 
are most strong in a polar crystal). The simplest and 
practical method of obtaining AtoC?) is to calculate the 
corresponding corrections to the single-particle energies. 
In the long- wavelength limit, upon comparing different 
physical processes and energy scales operative in the sys- 
tem, we can conceptually isolate, in specific properties, 
those effects arising from interactions between electrons 
and phonons from those caused by the other factors. This 
suggests Ato to be composed of three distinct elements, 
i.e. under the above assumptions one has in the low- 
doping regime: 



Ato(j) = (4iV)- 1 ^V(j,g) x A(j,q) x K{q) 
q 

where 
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A(k,q) = E 2 (k) - E^k + q) , g(qj) = gi 2 (qj) , and m 2 
denotes the band mass in the absence of the el-TO-ph 
interaction. The square of the TO phonon frequency 
fiq^ is represented as a sum of the unperturbed term 



Aw„ ; . We ex- 



n 2 

qj, which may be regarded 



the el-TO-ph interaction: 
pect the "normal" part w, 
as the bare TO phonon frequency with respect to the 
el-TO-ph interaction, to be free of softening anomalies. 
The proposed analytical decomposition, Eq. (j3|), allows 
us, firstly, to distinguish various physical factors con- 
tributed to AtoG?)) secondly, to treat them separately 
in order to understand how the interplay of these ele- 
ments affects the magnitude of Ato(j) ; and, thirdly, by 
using available experimental data, to obtain reliable up- 
per bounds for AtoO) ■ 

The first element V(j, q) is the main player of the de- 
composition since involves the effects of the hybridization 

of the electronic states by the TO phonon modes. It can 

1— 1 1 — 1 

be analyzed by using our previous results [7, 8] as follows. 
First of all, note that the quantity V(j, q) supplies a mea- 
sure of the zone-center TO phonon softening arising from 
the interband scatterings of the band electrons. That is, 
V(j, q) is proportional to a ratio of the electronic con- 
tribution I g(qj)| 2 /A(q), determined by the electron-ion 
potential and the band structure, and the lattice force 
constant MoW? associated with the bare TO vibrational 
mode. The momentum dependence of the quantity g(c{j) 
is seen from the first-principles analysis [8j: first, in a 
polar crystal, the maximum of g(qj) as a function of 
q is strictly related to the center point q = ; second, 
the function g(qj) is decreasing with increasing q. The 
last key aspect, which helps us to evaluate V(j, q) , is 
the delicate balance of long- and short-range forces pe- 
culiar to ferroelectrics (e.g., [l5|). By using the connec- 
tion between long-range forces and the el-TO-ph interac- 
tion |7|, |8j, one can relate the quantity g(qj) at q — >• to 
the relevant material parameters. As a result, in the long- 
wavelength limit, an upper-bound estimate for V(j, q) can 
be described by the expression: 



V(j,q) < 



S(j) 



3 > 



S(j) 



(7) 



Here, S(j) = S(j) /fi| and S(j) is a dipole oscillator 
strength associated with the jth zone-centre TO vibra- 
tional mode; is the core (electronic high-frequency) 
contribution to the dielectric function. In the derivation 
of Eq. ([7]) we supposed that contributions of the zone- 
centre TO phonons to the dielectric function e(o>) of the 
far-IR spectral r ang e are represented by a classical oscil- 
lator model 0, G3] : 



e(u,) 



S(j) 



^ fl 2 - 

3 U i° 



(8) 



It is seen that the parameter Tj always less than unity so 
that the inequality V(j, q) < 1 can be viewed as a reliable 
qualitative estimate. Note that, in a macroscopic sense, 
the difference between 1 and Tj can be considered as a 
measure of the polarity in such a manner that in polar 
materials, such as ferroelectrics, 1 
zero than in compounds of lower polar nature 



Tj is much closer to 
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Following 0, |||, we note that the second element 
A(j, q) serves as an adiabatic parameter which relates 
two different energy scales - one is the electronic scale 
determined by magnitude of an insulating gap E g ; the 
other is the conventional phonon frequency scale ftiOj . In 
view of the balance of long- and short-range forces as 
well as weak dispersion of the long-wavelength TO lat- 
tice vibrational modes, the product V(J, q) xA(j, q) can 
be taken to be momentum independent over the small q 
range. In the long- wavelength limit, this gives the esti- 
mate: A(j, Mljo/ (E g y/1 -t 3 - Mljo/Eg). 

The last element K,(q) associated with the derivates of 
A(k, q) measures at k — the effective curvature of the 
conductivity and valence bands. In contrast to V(j, q) 
and A{j,q), the role of the quantity lC(q) is not so im- 
portant for two reasons. First, it is not directly related to 
the lattice properties. Second, based on electronic struc- 
ture calculations in polar materials like ABO3 perovskite 
oxides, it is reasonable to assume that peculiarities of the 
relevant band structure close to the Fermi level as well 
as features of the Fermi surface topology do not critically 
influence the quantity lC(q). Therefore, for the sake of 
simplicity, assuming that the bands are parabolic around 
k = 0, we restrict here to the case K.(q)~ (1 + m 2 /mi) 
where mi corresponds to the bare effective mass in the 
valence band. (Note that in more general case the rele- 
vant estimate of /C may include the summation over q.) 
Of course, it is not guaranteed that such an approxima- 
tion of IC(q) would be almost the same for all the ABO3 
systems of interest. However, by above considerations, 
the assumption of the functional non- importance of !C(q) 
is quite evident in the present case. 

With these remarks in mind, the relatively simple ex- 
pression that can be used to find upper-bound estimates 
for AtoC/) i s obtained by gathering all three observed 
estimates together. The quantity Ato is then estimated 
in terms of observable parameters: 



AtoC?) 



< 



1 



1 



m 2 
mi 



\J 1 - Tj - Mljo/Eg E g 

(9) 

In Eq. ©, the adiabatic ratio Mljo/E g serves as a re- 
stricting factor which minimizes significantly the magni- 
tude of Ato(j)- This implies that there is no need to 
find more sharper estimates since the ratio of two pa- 
rameters, Ml jo and E g , is sufficiently small for systems 
under consideration. Thus, in order to ensure smallness 
of AtoO) f° r an Y polar compound of AB0 3 perovskite 
oxides, we need only the corresponding values of mate- 
rial parameters such as £00 , S(j), £ljo, and E g as well 
as the prefactor (1 + 7712/^1 ) (which can be deduced 
from electronic-band-structure calculations and/or from 
experimental data). Upper-bound estimates for the pa- 
rameter AtoO) calculated with the help of Eq. © for 
representative polar compounds belonging to typical per- 
ovskites are presented in Table 1. Note that these esti- 
mates can be referred to as first principles because Eq. © 
contains no adjustable parameters. 



TABLE I: Upper-bound estimates of AtoO) calculated for 
some typical representatives of polar perovskites and A IV B VI 
semiconductors. The estimates employ (1 + 1112/ mi) < 3 for 
the perovskites [J and (1 + m 2 /m 1 ) <2.2 for the A IV B VI 
compounds 18]. The rest of the material parameters were 
taken as in the previous work Q. 



Compound 


600 


E a (eV) 


3 


f!jo 


S{j) 


T j 


AtoO') 


Ato = 


SrTi0 3 


5.2 


3.3 


1 


88 


310 


0.977 


< 0.016 


< 0.02 








2 


176 


3.6 


0.333 


< 0.002 










3 


544 


1.6 


0.178 


< 0.003 




BaTiC-3 


5.2 


3.3 


1 


12 


1250 


0.994 


< 0.016 


< 0.02 








2 


182 


2.2 


0.234 


< 0.001 










3 


500 


0.8 


0.100 


< 0.002 




KNb0 3 


4.7 


3.1 


1 


96 


230 


0.972 


< 0.017 


< 0.03 








2 


198 


5.5 


0.451 


< 0.004 










3 


521 


3.1 


0.316 


< 0.006 




PbTe 


32.8 


0.31 




32 


632 


0.948 


< 0.031 




SnTe 


37.0 


0.35 




45 


330 


0.894 


< 0.025 




PbSe 


22.9 


0.28 




39 


289 


0.921 


< 0.033 




PbS 


17.2 


0.42 




67 


155 


0.890 


< 0.031 





Based on the comparative analysis of the obtained re- 
sults, we can draw several general conclusions and char- 
acterize the corresponding roles of the el-TO-ph inter- 
action as follows: 

(i) First of all, regarding effects associated with the 
el-TO-ph interaction, we stress a key distinction: On 
the one hand, it has been established that the large mag- 
nitude of the el-TO-ph interaction (~ 5 eV/A in the 
wide-gap compounds like ABO3 perovskite oxides and 
^0.6 eV/A in the A IV B VI narrow-gap semiconductors) 
is an important inherent property of polar crystals [7|, |8j . 
According to the vibronic theory (e.g., [3j), this provides 
a significant softening of the zone-center TO vibrations, 
and may lead to a structural instability of a crystal lat- 
tice. On the other hand, the contribution of the el-TO-ph 
interaction to A was found to be suppressed by the ratio 
Mljo/E g . For example, the data in Table 1 indicate that 
in the ABO3 materials the highest values of AtoC/) are 
in the range of 0.016-0.017 for the low-lying TO modes 
(j = 1). The high- lying modes have the smaller AtoC/) 
(j = 2, 3) because the relevant el-TO-ph interactions are 
smaller [7J. In addition to these numerical estimates, we 

have presented in Table 1 the corresponding results for 
a iv b vi 

narrow-gap semiconductors. It is easy to see that 
similar estimates are also obtained for these materials. 

(ii) The suppression of the contribution of the el-TO-ph 
interaction to the total electron-phonon coupling A is 
needed to be clarified. From a physical point of view, 
one may distinguish two factors why this is so. Firstly, 
note that one of the main features of the el-TO-ph inter- 
action is its direct connection with the relevant changes 
in the electron density distribution Q. The resulting 
shift of the electron density depends on interband transi- 
tions (hybridization) which exhibit an energy scale com- 
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pared with the magnitude of E g . In parallel, the carrier 
scattering effects associated with the long- wavelength TO 
phonon modes have the typical energy scale of the order 
of MljQ . Because the corresponding processes proceed in 
different ways, the smallness of the contribution of the 
el-TO-ph interaction to A can be explained by the great 
disproportion between these scales: MljQ<S^E g . Secondly, 
note that the el-TO-ph interaction directly depends on 
the electronic screening Q : the factor (eoo + 2) in Eq. ([7]) 
accounts for local-fields effects-induced partial screening 
of the bare el-TO-ph interaction. This implies that when 
the values of S(j) and are comparable, the local- fields 
effects are not to be ignored. 

(iii) The strength of the Frohlich interaction of electrons 
with the polar long-wavelength longitudinal optical (LO) 
vibrations can be measured by the coupling constant 
a. e ff represented as a sum of partial constants 
ot e ff = Si a(i) where the index i belongs to the rele- 
vant vibrational branches. A similar summation can also 
be introduced for the TO phonons as Ato = Sj Ato(j) ■ 
Since the effective contribution of the LO phonons to A 
is ~ a e ///6 (for weak coupling), the ratio ATo/(o ? e///6) 
gives then a qualitative comparison of the corresponding 
contributions associated with the TO and LO phonons. 
Using the polaron coupling constant a e //^2 reported 
in |S] and the data of Table 1, one can readily see that for 
SrTi03 the ratio Ato / ( a eff/6) < 0.06 . As it is apparent 
from this comparison, the polar LO phonons tend to pro- 
vide much stronger couplings of electrons than the long- 
wavelength TO phonons. This implies that the proper- 
ties of the conduction electrons are weakly affected by 
the the long-wavelength el-TO-ph interaction, no mat- 
ter how great the strength of this interaction is. 

(iv) In the context of superconductivity in SrTi03, it is of 
interest to note that, although the el-TO-ph interaction 
produces little influence on A , the relative contribution 
of the el-TO-ph interaction might be capable of playing 
some role in the doping phase diagram. (Here, we have 
followed the arguments presented within a two-band sce- 
nario of superconductivity fl9l. l20j|.) 

(v) Note that the derivation of Eq. (j9|) assumed an ad- 
equate choice of simplifications to recognize the relevant 
estimates for Ato ■ This allows us to test the applicability 
range of Eq. © to other systems. Let us demonstrate 
that Eq. © does also work in a Mott-Hubbard system 
iron monosilicide FeSi. Our goal here is with the help of 
the experimental data relating to FeSi to estimate values 
of Ato(j) for the four TO zone-center phonon modes and 
compare them with the corresponding values reported 
in [21|. Using the results of [22|], we choose the prefac- 
tor (1 + m2/mi)~ 2 within the qualitative picture of the 
electronic spectra. The calculated values of AtoO) are 
shown in Table 2 along with the values obtained in [2lj |. 
As expected, the agreement is fairly good. 

Note that the knowledge of material parameters, such 
as eoo, S(j), fljo, and E g , allows us to analyze the rele- 
vant interband hybridization effects in FeSi. In Ref. Q we 
have shown that this can be made by estimating the mag- 



TABLE II: Comparison of the values of Ato(j) for FeSi. 
The estimates were obtained from Eq. @ with the use of 
Eg = 0.14 eV (a direct gap value) £oo = 8 [23|], and 

experimental values of lattice TO vibration frequencies and 
strengths given in [2lj . 



TO mode (j) 


1 


2 


3 


4 




CljO (cm -1 ) 


198 


318 


338 


445 




S(j) 


8.0 


8.5 


2.8 


1.3 




Tj 


0.44 


0.46 


0.22 


0.11 


Ato = 


AtoCO 












This study 
Rcf. [21] 


0.068 
0.063 


0.142 
0.080 


0.056 
0.096 


0.041 
0.036 


0.307 
0.275 



TABLE III: Values of the el-TO-ph interaction constant g(j) 
and the transverse effective charge |Z*(J)| provided for each 
zone-center polar TO vibrational mode of FeSi. The estimates 
were made as described previously [?J with the use of the 
material parameters of Table 2. Two columns were added 
to allow comparison of results between different compounds. 
The entries in the columns labeled '(A)' and '(B)' refer to the 
estimates for FeO and PbTe Q, respectively. 



TO mode (j) 


1 


2 


3 


4 


(A) 


(B) 


g(j) ( eV/A) 


0.27 
1.63 


0.45 
2.70 


0.28 
1.65 


0.25 
1.47 


1.9 
2.19 


0.65 
8.37 



nitude of the effective constants g 2 (j) (representing the 
squared interband el-TO-ph interaction averaged over 
the given electronic states). The calculated values of 
g(j) for FeSi are given in Table 3. It is easy to see that 
the values of g(j) for both narrow-gap compounds ap- 
pear to be close to each other. In materials where the 
el-TO-ph interaction is operative, it can be verified by 
inspection of the magnitude of the zone-center polar TO 
vibrational mode effective charges |i?*(j)| 0, @|. This 
point is well illustrated in Table 3, which presents the es- 
timated values of the transverse effective charges, |Z*(j)|, 
both for wide-gap (FeO) and narrow-gap (FeSi, PbTe) 
compounds. Thus, the estimates presented in Table 3 
suggest that the interband channels of electronic scatter- 
ings associated with the zone-center TO phonons play a 
certain role in the lattice-dynamical properties of FeSi. 



IV. CONCLUSION 

In the present paper, motivated by the need to extend 
the results of the previous works @, Q , we have analyzed 
the influence of the interaction of electrons with infrared- 
active long-wavelength TO lattice vibrations on the to- 
tal electron-phonon coupling constant A . To examine in 
detail in what extent el-TO-ph interactions are impor- 
tant for the total strength of electron-phonon coupling we 
considered the generic two-band model involving the in- 
terband hybridization by the TO vibrational modes. The 
main result of the analysis points out the marginal role 
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played by the el-TO-ph interaction on electron pairing in 
doped polar insulating systems: the presented analytical 
and numerical arguments show that the partial contribu- 
tion from the el-TO-ph interaction to A is much smaller 
than the relevant contribution from the interaction with 
the polar long-wavelength LO phonons. The correspond- 
ing explanation is that the important physical properties 
of the el-TO-ph interaction are described by the matrix 
elements which mix the electron states of the valence 
and conduction bands separated by relatively large en- 
ergy gap. For given intra- and interband interactions 
coexisting in a doped insulator, due to distinct energy 
scales, the amount of the contribution of the el-TO-ph 
interaction is ruled by the adiabatic ratio HQo / E g . In 
most systems of interest, the quantities, E g and hQjo, 
are of a different order of magnitude so that the real im- 
pact of the el-TO-ph interaction on A turns out to be 
strongly suppressed. This conclusion has a clear physical 
meaning: upon completion of the intraband scatterings, a 



conduction electron becomes already nearly "undressed" . 
Our studies suggest that the el-TO-ph interaction (via 
valence electrons) is of great concern for the lattice and 
polarizability properties, the changes in the electron den- 
sity distribution, and the renormalization (soft behavior) 
of the zone-center TO vibrational modes. At the same 
time, in the context of superconductivity in doped po- 
lar insulating compounds, the el-TO-ph interaction due 
to a difference in energy scales of underlying interactions 
tends to contribute little to the total effective coupling 
constant. 
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